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1. Introduction 

The achievement of Bose-Einstein Condensation (BEC) in dilute gases has been quickly 
followed by the demonstration of atomic outcouplers |U [21 El H], allowing atoms to 
be extracted coherently, all in the same wavefunction, forming a so-called atom laser. 
Continuous outcoupling [31 H] leads to quasicontinuous atom lasers, and among them 
gravity compensated guided atom lasers (GAL) [3 El [7] are promising tools to study 
quantum transmission through obstacles [El El [101 [HI [El [13], as they provide a large 
and constant de Broglie wavelength along the propagation. Tuning the atom laser flux 
should allow one to investigate both linear and nonlinear phenomena, and even reach 
strongly correlated ones [8]. However quasicontinuous emission imposes to operate 
in the so-called weak coupling regime (see e.g. [H] and references therein) such that 
corresponding flux limitations may restrict the envisaged applications. Thus, this work 
aims to characterize in detail the flux properties of the horizontally, radiofrequency (rf ) 
outcoupled, GAL demonstrated in [5], in order to identify the limits of the weak coupling 
regime. 

The weak coupling limits are intimately linked to the irreversibility of the coupling 
process; the coupling strength has to be weak enough for the extracted atoms not to be 
coupled back into the trapped BEC [151 EH EU] ■ In particular, the emergence of bound 
states as the coupling strength is increased constitutes a major limitation [HI [19] and can 
even result in the shut-down of the atom laser [201 [2~Tj . In this respect, the crucial role 
played by the underlying force in the coupling process (arising either from the gravity 
acceleration in the case of free space rf outcoupling or from the two-photon momentum 
kick in the case of Raman outcoupling) has been stressed [13] ■ This force pushes out 
the atoms and prevents them to be recaptured before they have left the coupling zone. 

For the horizontally GAL studied here, gravity is cancelled, and the extracting 
force is only provided by the repulsive atomic interactions with the remaining trapped 
BEC. This has important consequences on the flux properties. First, the outcoupling 
spectrum, i.e. the dependence of the output flux versus the rf frequency, is profoundly 
modified. From a calculation based on the Fermi golden rule, we indeed show that 
it exhibits a double structure, with a large pedestal and a very narrow peak around 
its maximum. This is in stark contrast to rf-outcoupled free falling atom lasers [3]. In 
practice, it results in a strong sensitivity to residual magnetic fluctuations, which should 
be suppressed as much as possible. Second, the horizontally GAL is subjected to severe 
flux limitations, as expected from the weakness of the extracting force. In particular, we 
explicit here the link between such limitations, arising from the onset of bound states 
in the process as recognized in other schemes [201 EH EH], and the so-called "quantum 
pressure" associated with the confinement of the atom laser wave function around the 
outcoupling zone. Altogether, the maximal flux achievable is typically one order of 
magnitude less than for the gravity driven atom lasers. However these limitations do 
not impair future applications for studying quantum transport phenomena, as we finally 
discuss. 
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Figure 1. Generation of an atom laser by rf outcoupling from a BEC created in an 
hybrid opto-magnetic trap. The optical guide is horizontal such that the atoms do 
not undergo any acceleration after leaving the BEC zone. Note that an atom laser is 
emitted on both side. For clarity, only one is shown here. 



Both aspects (the outcoupling spectrum and the weak coupling limits) are 
investigated using a simplified quasi- ID model, whose predictions are compared, with 
an overall fair agreement, to the experiments. The paper is organized as follows. We 
present in section |2]the scheme of the GAL demonstrated in [5] and the quasi-lD model 
we use. Then we study in section [3] the outcoupling spectrum and detail in section @] 
the weak coupling limits. Finally we conclude and discuss some prospects in section |5j 

2. The rf guided atom laser: principle and modelling 

2.1. Geometry and outcoupling scheme 

As depicted in figure (TJ a BEC containing iVo 87 Rb atoms is directly created inside an 
optomagnetic trap in the \F,m F ) = 1 1 , — 1 ) magnetic state. This cigar-shaped trap 
is the combination of an horizontal optical guide, which generates a strong cylindrical 
transverse confinement (with a radial frequency u>±), and a quadrupole magnetic field, 
which provides a weak longitudinal confinement (of frequency u z <C u±) along the z 
axis. The chemical potential of the BEC verifies /^bec 3> ^>±, such that the BEC is well 
in the 3D Thomas- Fermi (TF) regime [22]. Its density nsi)(rj_,z) (r*j_ being the radial 
coordinate) has then an inverted parabola shape of radii R± iZ = (2/iBEC / ' mw i i2 ) 1/ ' 2 > with 
m the atomic mass. 

The atom laser is generated using an rf magnetic field to coherently spin-flip a small 
fraction of the atoms from the source BEC to the |1, 0) state which is insensitive (at first 
order) to the magnetic field. While the outcoupled atoms are still radially confined by 
the optical guide, they can move along the longitudinal axis z. More precisely, they feel 
the sum potential V g (r±, z) + gn 3 n(r^, z). Here V g takes into account the optical guide 
confinement together with the second order Zeeman quadratic effect, and the second 
term corresponds to the mean-field interaction with the remaining trapped BEC. Here 
g = Airah 2 jm is the collisional parameter and a the scattering length. 

This interaction provides the necessary force to extract the atoms from the coupling 
zone. The atoms accelerate down the "bump-shaped" mean-field potential until they 
leave the BEC region and then propagate at constant velocity v, i.e. constant de Broglie 
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wavelength AdB = h/mv [5]. 

2.2. Hypothesis and quasi- ID model 

In the following, we use a simplified model to describe the GAL. It relies on the 
assumptions listed below: 

(i) The influence of the |1, 1) magnetic state is negligible and we only consider a two- 
level system (see section I3T4")) . 

(ii) The BEC wave function evolves adiabatically during the outcoupling process. Apart 
from section 14.31 we also neglect the depletion of the BEC. 

(iii) The collisional parameter g is independent of the magnetic substate, which is 
approximatively the case for 87 Rb [23J. 

(iv) The intra-laser interactions are negligible in the outcoupling process. Note that this 
does not prevent the interactions to play a significant role along the propagation in 
the guide (see section [SJ). 

(v) The atom laser does not experience any residual force along the guide axis z: 
Vg(r_\_, z) = mui 2 _r 2 _/2. This can be achieved by using the weak longitudinal optical 
guide confinement to compensate precisely the second-order Zeeman effect [5]. 

Besides, we adapt the usual quasi-lD approximation (see e.g. [24J) to our 
configuration where two magnetic substates, the atom laser and the BEC, are involved. 
If one consider only one single state, the quasi-lD approximation supposes that the 
transverse wavefunction adjusts adiabatically along the guide axis: it sticks to the 
ground state of the local transverse potential, which depends on the total mean-field 
interaction, i.e. including both the inter- and intra-states contributions. In our case, the 
hypothesis of g being independent of the magnetic substates is of particular importance: 
it implies that the atom laser and the BEC experience the same transverse potential. 
Their transverse wave function should thus be identical, and the same as if all atoms 
were in the same state. Consistently, we assume a perfect transverse mode matching 
and write the total wave function i/?(r, t) in the separate form 







(j) B (z,t) 









^±{r±,z), (1) 



with the normalization condition J \if)\ 2 drj_dz = 1. Here the subscript B and L refers 
to the BEC and the atom laser respectively. Neglecting the intra-laser interactions, ip± 
evolves continuously from the transverse inverted parabola shape of the BEC to the 
gaussian fundamental mode of the guide. 



As detailed in |Appendix A[ the outcoupling process is then reduced to a one- 
dimensional problem. In particular, the effective ID mean-field potential felt by the 
atom laser follows the BEC profile (see figure [2]) and reads 

2 I I 

Mi (*) = &"± + Abec (i-t)e(i- |i) (2) 
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the origin of the z-axis being at the center of the BEC. Here G is the Heaviside step 
function, /Ibec = /Ubec — is the chemical potential corrected from the transverse 
fundamental mode energy and R z = (2p,BEc/ m ^z) 1 ^ 2 - The latter corresponds precisely 
to the BEC edge in the quasi-lD approximation. In the following, we omit the constant 
in (j2J) and use the implicit energy shift E — >■ E — hu±. 

3. Outcoupling spectrum: flux vs rf frequency 

We study now the dependence of the flux versus the rf frequency, i.e. the outcoupling 
spectrum. In particular, we show that it exhibits a very peaked shape. Note that we 
assume throughout this section that the atom laser operates within the weak coupling 
regime, which is defined below. The associated conditions will be detailed in the next 
section. 

3. 1 . Theoretical framework 

As sketched in figure El the rf outcoupling process can be seen as the coupling of an 
initial state (BEC) to the continuum formed by the stationary states in the potential 
Vji(^). We consider here the weak coupling regime, which is defined by the regime 
where the coupling process is irreversible: the extracted atoms leave the coupling zone 
irretrievably, without any coherence with the remaining trapped BEC [T6 | fT5 | fT7]. It is 
also commonly referred as the Born-Markov regime [25]. In this situation, the emission 
is quasicontinuous and the BEC state decays exponentially, the outcoupling rate Y (or 
equivalently the flux J 7 = N Q T) being given by the Fermi golden rule [17] 




(3) 




BEC 



Z E R Z 



Figure 2. Effective ID energy diagram for the outcoupling process. V\\(z) 
[equation ([2])] is the potential seen by the atom laser state, /ibec = A*bec — fUo± is 
the corrected chemical potential and R z the BEC edge. 5u T f is the frequency detuning 
from a coupling at the center of the condensate and ze corresponds to the classical 
turning point given by equation ([5]) (only the z > resonance is shown). 
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In this expression, t](E) is the ID density of states and fi r f quantifies the rf coupling 



strength (see Appendix A). </>b corresponds to the longitudinal BEC wave function 
normalized to unity. It is well approximated by |0b(z)| ~ (15/16i? 2 ) 1//2 (l — z 2 /R 2 z )Q(l — 
z/R z ) except around the BEC edge (z = R z ), where it falls down to zero. Lastly, 4>i : e 
is the atom laser stationary wave function, also normalized to unity, whose energy E is 
selected by the rf frequency u r f/2n via the resonance condition 

E = E BEC - fiwrf . (4) 

The calculation of T(E) according to (j3J) reduces then to the determination of 0l,e (see 



Appendix B ). However, in most semiclassical approach is sufficient to obtain 



an accurate approximation |17j . 
3.2. Semiclassical approach 

Following the so-called Franck-Condon principle, we consider here that the outcoupling 
is spatially localized at the points defined by energy and momentum conservation during 
the coupling, i.e., the classical turning points [26J. As both the initial kinetic energy of 
the condensed atoms and the momentum transfered by the rf field are negligible, the 
turning points are defined by V(ze) = E. It gives the two symmetrical positions 



2h5u ri 

V mUJ z 

where 5u r f = u r f — (Ebec — /Wbec)/^ is the frequency detuning from a coupling at 
the center of the condensate (see figure [2D- Formally, this approach is equivalent to 
approximate the atom laser wave function fa^z) by a Dirac distribution around Ze, 
such that equation (j3J) transforms into [27] 

T sc (E) = ^Si j dz\Mz)\ 2 S{V\\{z)-E). (6) 

Considering the emission at one side, it finally reduces to the same expression as in [TF] 

the gravity acceleration being replaced here by the local acceleration g e g = lo\ze of the 
potential V\\ at the position ze- Consistently it vanishes for a coupling at the BEC edge 
(ze = R z or equivalently 5cu T f = /2bec/^)) the corrected chemical potential being the 
outcoupling spectrum width in the semiclassical approach (see figure [2]). 

This approach is well justified when one outcouples at the BEC side. There the 
acceleration is large enough and the first lobe of the atom laser wave function (which can 
be approximated by an Airy function) is narrow compared to the size of the BEC (see 



Appendix Bp . However it fails near the top of the BEC (g e g — > 0), where the coupling 
rate should precisely reach its maximum. This region is thus of particular interest and 
calls for a more complete description. 



Quasicontinuous horizontally guided atom laser: coupling spectrum and flux limits 7 

8 



6 



s!S 4 
U 

2 


0.0 0.2 0.4 0.6 0.8 1.0 

hSa> rf I p BEC 

Figure 3. Outcoupling rate versus the normalized detuning Mw r f//iBEC- Straight 
(red) line: outcoupling rate T given by the Fermi golden rule [equation ||3J)]. Dashed 
(blue) line: semiclassical approximation r sc [equation ([7|l]. Inset: zoom on the 
vicinity of the zero detuning (the normalization is changed to uj z ). The parameters 
(Abec = 2.5 kHz, oj z /2tt =27 Hz and fi r f/27r = 15 Hz) are those of section EHtl 




3.3. Outcoupling spectrum 

The determination of the stationary atom laser wave function 0l,e is detailed in 



the |Appendix B[ Then, using equation ([3]), we calculate the outcoupling rate. The 
corresponding spectrum T(Su r {) is plotted in figure El 

As expected the outcoupling spectrum extends over a typical width /ibec/^ and 
agrees well with the semiclassical expression r sc when coupling at the BEC side 
(5u) r f ^> u z ). De facto, its validity range is very broad and a significative deviation only 
appears as 5u T { ~ u z . Such behaviour can be qualitatively understood noticing that the 
classical turning point ze approaches here the characteristic length a z = (h/muz) 1 ^ 2 
of the inverted harmonic potential [see equation (JSJ)]. Thus the curvature dominates 
and leads to the saturation of the first lobe of the atom laser wave function L e (see 



Appendix B[ ). We thus expect the outcoupling rate to reach a maximum around this 
position. Finally, it should fall down rapidly to zero for 5oj t { < 0, as it corresponds 
to the classically forbidden region. Such analysis is refined by a zoom on this region 
(see the inset of figure [3]). In particular, it appears that the maximum value is reached 
around 5u r { ~ oj z /A and 

max(r)~r sc (^) = ^n CTz . (8) 

ow rf V 4 / UJ Z 

Here IL^ = 7icr z \(f) B (0)\ 2 /\/2 corresponds to the proportion of atoms contained in the 
region of spatial width -ko~ z /V2 0. 



X Consistently, we recover the expression for a coupling to a continuum of finite width hui z , the Rabi 
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Altogether, the outcoupling spectrum exhibits a double structure. The large one, 
of width given by the chemical potential /iBEC; corresponds to a coupling at the BEC 
side. The very narrow one, of width u z , is associated to the sharp increased of T around 
the BEC center. For a BEC in the 3D Thomas- Fermi regime (/Ubec ^ ^J- ^ ^z), it 
results in a very peaked outcoupling spectrum, in contrast with the common behaviour 
of free falling atom lasers (see e.g. [TT]). 

Such peaked shape, whose narrow width is typically around a few tens of Hertz, 
has important consequences for the use of such GAL. Indeed weak residual magnetic 
fluctuations are hardly avoidable in practice. First those will slightly broaden the 
spectrum such that the narrow peak should be hard to observe experimentally. Most 
importantly, those may result in strong density fluctuations when coupling near the 
BEC center. This indicates that an extreme attention should be paid when operating 
in this region. As we shall see in section HJ this statement is even more reinforced as 
this position corresponds to the most severe weak coupling conditions. 

3.4- Experimental investigation 

We present now our experimental investigation on the outcoupling spectrum. The 
optical guide is created by a Nd:YAG laser (A=1064 nm) focused on a waist of 
30 /im, resulting in a radial trapping frequency u±/2tt = 330 Hz. The weak magnetic 
longitudinal confinement is set to u z /27r = 27 Hz. The initial BEC contains N = 
1.1 x 10 5 atoms, corresponding to /Ubec / h = 2.5 kHz. 

The bias field of our Ioffe-Pritchard magnetic trap is chosen around B ~ 7 G 
(see [29] for further details). For such bias, the Zeeman quadratic effect enables us to 
perform selectively the rf coupling to the atom laser state, without populating the |1, 1) 
state [30]. We calibrate the coupling strength by using very short outcoupling pulses 
(t r f "C h/^BEc)- There the whole BEC oscillates precisely at the pulsation f2 r f, which is 
determined within a few percent. 

In contrast, a quasicontinuous GAL is extracted from the BEC by applying a 
long and weak rf pulse. In the following, the coupling time is set to t r { = 50 ms. 
It should allow us to resolve, in principle, the peaked shape of the outcoupling spectrum 
(trf 1/uj z ^> ft//iBEc)- An experimental spectrum is shown in figure H] for a 
coupling strength Q t {/2tt = 15.1 Hz. The extracted atom number for one side, N-^, is 
measured by fluorescence imaging and used to calculate the outcoupling rate according 
to T = N L /N t T f, the depletion of the BEC being negligible for our parameters. 

For a coupling at the BEC side, we find a fair agreement with the predictions 
(without any free parameters). For a coupling near the BEC edge, we observe a slight 
deviation that can be traced to the presence of residual thermal atoms. This broadens 
the spectrum, leading to a total width (~3 kHz) slightly above the ideal value given by 
the corrected chemical potential /Ubec/^- 

Looking at the coupling near the zero detuning, the situation is balanced. On 



oscillation frequency being fl r [ 
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Figure 4. Experimental outcoupling spectrum T(Su I {) for the coupling strength 
fi r f / 27r = 15.1 ± 0.5 Hz. Left vertical axis: outcoupled atom number (for one 
side). Red circles: experimental results. Straight line: theoretical predictions given by 
equation [3J The total atom number is A*o = 1.1 x 10 5 and /^bec/^ =2.5 kHz. The 
horizontal error bars (200 Hz rms) corresponds to the uncertainty on the detuning 
caused by the shot-to-shot magnetic fluctuations. Note that the detuning has been 
corrected from the long term bias drift estimated around 0.4 mG/h (300 Hz/h). 

the one hand, we recover a peak-shaped structure associated with a sharp variation of 
the flux. It shows that the rapid magnetic fluctuations have been fairly suppressed. 
In particular, one can infer an upper bound around 70 fiG rms (50 Hz converted in 
frequency unit). This was achieved by using a magnetic shielding and actively stabilizing 
the current Iq that creates the bias field (AI/I ~ 10 ppm rms). 

On the other hand, the spectrum is "scrambled" in this region by the remaining 
shot-to-shot fluctuations. Those can be estimated around 250 yumG rms (200 Hz) and 
prevent at the moment a detailed investigation of the narrow peaked structure. De 
facto, the observed maximum outcoupling rate lies around T ~ 1.5 s" 1 , still far from the 
expected maximum (see figure [3]). As we shall see in the next section, let us stress that 
this difference may very unlikely be attributed to saturation effects beyond the weak 
coupling regime. 

4. Flux limitations for the rf guided atom laser 

In this section, we explicit the weak coupling limits for the GAL. We show that they 
originates from the emergence of bound states in the coupling process and we especially 
emphasize their close link with the "quantum pressure" for rf coupling schemes. 
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4-1- Weak coupling conditions 

As already mentioned, the weak coupling regime is referred as the situation where the 
coupling strength is low enough for the dynamics to be irreversible, i.e. well described 
by the Born Markov approximation. This imposes several conditions associated with 
different physical phenomena. 

Building upon the modeling of the outcoupling process as the coupling to a 
continuum, one can first derive a necessary condition: the memory time t m of the 
continuum has to be short enough for the extracted atoms not to be coupled back to 
the trapped state, i.e. Tt m <C 1 [161 US]- This finite time t m being linked to the finite 
effective width of the continuum, t m ~ h/A, the condition can be equivalently written 
in the canonical form r < A [28] [T7] . From the double structure of the outcoupling 
spectrum discussed in section 13.31 we expect this limit to evolve from r < w 2 for an 
outcoupling near the BEC center to T <C /^becM a t the BEC side. 

However, a more strict condition arises from the emergence of bound states in the 
rf dressed basis, as first mentioned by Stenholm et al. [311 E2] in the context of laser 
induced molecular excitations. If the coupling strength is increased above a certain 
critical value, part of the atoms remains indeed trapped in a coherent superposition of 
the atom laser and BEC states, and do not escape the coupling region. Such phenomenon 
leads to a strong modification of the outcoupling dynamics [HI [19] and ultimately results 
in the shut-down of the atom laser, as observed experimentally [201 12H E]. 

Following the analysis of Debs et al. [H], the critical coupling strength VL c rf coincides 
with the adiabatic following condition in the rf dressed potentials. These potentials are 
shown in figure [5] for an outcoupling at the BEC side. In this case, the local acceleration 
<7eff = uj^ze is finite at the anti-crossing position [i.e. the classical turning point Ze given 
by equation fl5])] and the results of Ref. [33] apply. From the transfer probability between 




Figure 5. Energy level crossing in the dressed state basis at given rf detuning 5aj r [ and 
coupling strength i7 r f. Dashed lines: bare potentials for the BEC (|1,— 1)) and atom 
laser (|1,0)) states. Solid lines: dressed state potentials. The crossing position ze 
corresponds to the classical turning point and g e s = ui^ze the local acceleration. The 
typical anti-crossing extension is given by Azn rf = ^rf/ m 9cff- Red line: stationary 
atom laser wave function 0l,e- £z e is the typical extension of the first lobe (see text). 
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the two dressed states [§[ one can therefore infer the weak coupling limit 



for 5uo T i ^> Co> 2 . Such expression enlightens particularly well the importance of having 
a large local acceleration, that pushes out the extracted atoms and prevents their 
recapture, to achieve a large critical coupling, i.e. a large flux [TJ]. 

Interestingly, one finds that the anti-crossing extension (see inset of figure [5]) 
corresponding to the critical coupling strength Azqc = hVL c vi /mg e R coincides precisely 
to the typical Airy lobe size C ZE of the stationary atom laser wavefunction 0l,e : 
Az n «, =(h 2 /2m 2 g cS ) 1 / 3 = C ZE . 

This suggests an alternative picture to recover the limit given by equation ([9]). As 
a matter of fact, the dressed state basis diagonalizes the bare potentials and coupling at 
each position, i.e. the complete Hamiltonian, except for the kinetic energy. Formally, it 
means that the kinetic energy term drives the coupling between the two dressed states. 
The decay from one state to the other is then governed by its relative importance 
compared to the coupling strength (see e.g. [33]). The adiabatic following condition 
can thus be reformulated as E}.(ze) <C hVL c {t , where Ek(zE) is the kinetic energy at 
the crossing position. In the case of a rf outcoupler, no momentum kick is imparted 
to the atoms and, in a semiclassical picture, they leave the outcoupling zone without 
initial velocity. However E^ze) ~ h 2 /2m£ 2 E remains finite and is given by the so- 
called quantum pressure associated to initial confinement. This leads directly to the 
criterion ([9]) (apart from a scaling prefactor (2/7r) 2//3 close to unity). 

For an outcoupling near the BEC center, the local acceleration vanishes and the 
validity of (02) breaks down. However, the above approach based on the quantum 
pressure can be directly extended in this region, the Airy lobe size being replaced by 
the harmonic oscillator size o~ z . The critical coupling writes then 

°& ~ (2^75 ( 10 > 

for &Lo x i < u z (for consistency, the scaling prefactor found above has been included). 
Note that it corresponds to the limit found in [181 IT9] . considering non- interacting 
atoms coupled into free space. 

Altogether, we bridge the gap between the two domains using the phenomenological 
expression 



(2vr 2 ) 1 /3 



5u 



^ f (^ rf ) = -^— i + 64(^) (11) 



UJ Z 



1/6 



over the whole coupling spectrum [see figure [6] (a)] 



4-2. Critical outcoupling spectrum 



From the limit given by equation ( ITTjl . we deduce the critical outcoupling rate T c = 
r sc (^rf > ^rf) for each detuning. Note that we use here the semiclassical expression (J7J) 

§ The transfer probablity P - exp - (n ri /n c ri ) 3 / 2 of Ref. [33] can be seen as the extrapolation of the 
Landau-Zener formula for a particle starting at zero velocity at the crossing position. 
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Figure 6. (a) Evolution of the critical coupling strength f^f with the normalized 
outcoupling detuning Mw r f / /Ibec- Straight and dashed lines: f2£f given by 
equations, (jlip and ([9]) respectively. Inset: zoom around the zero detuning, with 
a normalization changed to u> z . (b) Critical outcoupling spectrum r c (5w r f) = 
r sc (f^ f , 6u) T { ) for the detuning range Suj t { £ [u> z /4, //bec/S-]- Tmax is the maximal 
coupling rate and Suj™ ax the associated detuning. The parameters are the ones of 
section [X^l uj z /2tt = 27 Hz and /}bec/^=2.5 kHz. 



and restrain the detuning range to 5u T f G [w*/4, /2bec/^] where this approximation is 
accurate. The resulting critical outcoupling spectrum is plotted on figure[6] (b). Here the 
very peaked structure has disappeared: the sharp decrease of r sc (5o; r f) is indeed partly 
compensated by the concomitant increase of Q^(Scu T {). In particular the evolution of T c 
is rather smooth around its maximum T max , i.e. the overall maximum outcoupling rate 
achievable within the weak coupling regime. Such maximum reads 



the optimal outcoupling position <5w™ ax ~ 0.08 /2bec/^ being shifted towards the BEC 
edge 1{J. 

The above upper bound (fT2]) is of particular importance as it allows us to explicitly 
define the operating range of the atom laser. It imposes a strict limitation on the 
outcoupling rate, having in particular r max <C uj z . Considering the experimental 
parameters of the previous section, one finds r max ~ 8 s _1 , i.e. a maximum flux 
Jnax = Tmax ~ 10 6 at.s -1 . This is typically one or two orders of magnitude less 
than with "common" schemes (see e.g. [33]). As already said, such low flux operation 
was expected for the GAL. It can be traced to the very weak extracting force (only due 
to the interaction with the BEC) experienced by the atoms. 

Finally, we can estimate the critical coupling strength for a detuning 5u T { ~ w 2 /4 
corresponding to the maximum outcoupling rate in the weak coupling regime (see 
figure [3]). As shown in the inset of figure [6], we find ~13 Hz around this position, 
leading to T c ~ 6 s _1 . The coupling strength used to calculate the outcoupling 
spectrum presented in figure [3] being slightly above such critical value (O r f = 15.1 Hz), 

|| Note that r(<5w r f) being less sensitive to residual fluctuations in this region (see section l373|) . the 
corresponding detuning <5w™ ax is also optimal from a practical point of view. 




0.22 



(yU B Ec/^) 1/3 ' 



(12) 
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effects beyond the weak coupling regime should a priori be taken into account in the 
predictions. However those may unlikely explain the experimental observations shown 
in figure HI the apparent saturation being most certainly due the residual shot-to-shot 
magnetic fluctuations. 



4-3. Experimental investigation 

Here we investigate experimentally the flux limitations by monitoring the extracted atom 
number Al (still for one side) as the coupling strength f2 r f is progressively increased. 
The trap frequencies are the same as in section 13.41 and the initial atom number is now 
iVo = 2 x 10 5 atoms, leading to pBEc/h = 3.3 kHz. The outcoupling time is decreased to 
£rf = 20 ms in order to keep the BEC depletion small enough for large outcoupling rate. 
The results are shown in figure [7] for three different detunings (<5a; r f/27r = 0.6, 1.3 and 
2.3 kHz). Note that they correspond to outcoupling positions at the BEC side, where 
the extracted atom number is less sensitive to the magnetic shot-to-shot fluctuations. 

As expected jV L agrees with the semiclassical predictions for weak coupling 
strengths (see inset of figure [7J. However, a significative deviation appears when O r f 
is increased and starts to saturate. If a part of the deviation can be imparted to 
the depletion of the BEC, it cannot fully explain the observations. This is confirmed 
by comparing the experimental results with the solution El of dNL/dt = F sc (Nb)Nb, T sc 
being still given by equation ([7]) and Ab — — 2Al being the BEC atom number. 
In contrast, the numerical resolution of the coupled quasi-lD Gross-Pitaevskii (GP) 
equations [equations. (IA.5j) and flA.6j) ] fits well the experimental data, especially for 



the detuning <5z/ rf = 2.25 kHz (red curve). The saturation at large Al for the detunings 
5u r{ = 0.6 and 1.29 kHz (green and blue curves) will be discussed below. 

The disagreement between the experiments and our calculations based on the Fermi 
golden rule can be interpreted clS du SI gnature of the breakdown of the weak coupling 
regime. As a matter of fact, the onset of such deviation is well rendered, without any 
adjustable parameters, by the limit defined by equation ( [Til . Thus, these observations 
comfort the analysis presented in section 14.11 as previously stressed for other schemes 
[20l I2T1 [14] , the predominant flux limitation likely originates from the presence of bound 
states in the rf dressed potentials. As already said, equation ( TTTj) finally allows us to 
properly define the quasicontinuous operating range of the guided atom laser, which is 
symbolized by the gray area in figure [7J In particular, it gives a maximal Al around 
3 x 10 4 atoms for the parameters considered in this section, corresponding to a maximal 

flux Jmax ~ 1.5 X 10 6 at.S _1 . 

Before concluding, one should stress that these results support a posteriori the 
approximations made to model the GAL and to calculate the outcoupling rate. In 
particular we found no signatures (on the flux) of any BEC excitations in the weak 

% Using either the Fermi golden rule ([3]) or the semiclassical expression ([7]) gives very negligible 
differences for our parameters. 

+ An analytical solution can be derived as detailed in |34j . 
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Q rf /2;r (Hz) 

Figure 7. Experimental investigation of the weak coupling limits: the extracted 
atom number for one side vs the coupling strength f2 r f for different detunings. 
A (green): 5u> vl /2ir = 0.6 ±0.1 kHz. o (blue): <5w rf /27r = 1.3 ± 0.1 kHz. o (red): 
(5w r f/27r = 2.3 ± 0.1 kHz. The outcoupling time is t T { = 20 ms. Dash-dotted lines: 
semiclassical prediction Ni, = r sc (&j r f) i r f . Dotted lines : the same with the depiction 
of the BEC being taken into account (see text). Thick lines: numerical resolution of 
the coupled quasi-lD GP equations. The grey area corresponds to the weak coupling 
regime defined by equation (|11|) . Its border (thin line) is given by the parametric curve 
[iV£(&j r f), f2£ f (&j r f)], Nl = r c (5oj T i)t T f being the critical outcoupled atom number. 
The inset is a zoom, in log-log scale, on low rf coupling strength. 



coupling regime defined by ffTTj) (grey area in figure [7]). The same applies for 
the intra-laser interactions. For the latter, their influence is characterized by the 
dimensionless parameter an^, where is the atom laser linear density. In the guide, 
it can be estimated using n L ~ J- /v. Considering the typical propagation velocity 
v ~ (2/iBEc/' m ) 1 ^ 2 ~ 5 mm.s -1 together with the maximal flux J-" max , one finds 
(cmijmax ~ 1-5, i.e. slightly above the limit of the ID mean-field regime defined by 



an L ^ 1 (see Appendix A ). Increasing the flux further, the intra-laser interactions 
will profoundly modifies the atom laser transverse wavefunction and so the transverse 
dynamic. Then our quasi-lD model may fail and it could explain the differences 
observed between the numerical simulations and the experiments at large outcoupled 
atom number (green and blue curves in figure [7]). 

5. Conclusion and discussion 

We characterized in this paper the flux properties of the guided atom laser. In particular 
we showed that a quasicontinuous emission entails severe flux limitations (J 7 < 10 6 
at.s -1 ) compared to other schemes. However, it does not a priori alter the very good 
promises for studying quantum transport phenomena, both in the linear and nonlinear 
regime. The inherent weak flux is indeed balanced by the very slow velocity propagation 
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(and thus large de Broglie wavelength) offered by the rf outcoupling. It results in not so 
dilute atomic beams, typically within the ID mean-field regime (ani, < 1) as discussed 
above. In addition to single particle quantum effects (for instance quantum tunnelling 
or Anderson localization) a rich variety of nonlinear phenomena can be observed, as for 
instance the atomic blockade or a soliton generation past a finite lattice [10] . 

As a matter of fact, nonlinear effects may show up for very dilute atomic beams, 
as well illustrated by the transmission through a ID disorder slab (see e.g. [T2] ). There 
the interaction-induced localization-delocalization crossover should be observed for the 
critical sound velocity c = (2hjj±an^/m) 1 ^ 2 ~ v/7, i.e. deep in the supersonic regime. 
Taking v = 2 mm.s -1 as in [35], it corresponds to an atomic density an^ ~ 0.05 (in 
dimensionless unit), well within the limits given in this paper. 

In future, the scheme described here could be extended to other configuration. 
Along these lines, one may implement Raman outcoupling techniques or head towards 
"all optical" techniques, similar to the one described in [6] but using Bragg transitions 
as an output coupler [36J. Besides, tuning the scattering length through Feshbach 
resonances (using a proper atomic species) would certainly offer new possibilities (see 



To conclude, let us stress that the complete characterization of the rf outcoupled 
guided atom laser requires a detailed study of its longitudinal coherence, i.e. its energy 
linewidth. This will be the subject of a future work. 

Acknowledgments 

We thank J.-F. Riou for fruitful discussions. This research was supported by CNRS, 
Direction Generale de l'Armement, ANR-08-blan-00 16-01, IXSEA, EuroQuasar program 
of the EU. The Laboratoire Charles Fabry de l'lnstitut d'Optique is a member of IFRAF. 

Appendix A. Quasi-ID model 

We detail here the quasi-lD model for the guided atom laser. First we derive the quasi- 
1D coupled GP equations for the longitudinal wave functions 0b (z, t) and (fih(z, t). 
Second we present the approximated form resulting from the assumptions listed in 
section 12.21 

Appendix A. 1. Quasi-ID coupled Gross-Pitaevskii equations 

The rf magnetic field B r { cos w r f t (along a direction perpendicular to the biais field B ) 
couples the BEC and atom laser states. In the rotating wave approximation, it results 
into the "mean-field" Hamiltonian 



e.g. [37]). 



2 



(Al) 




A + v g + iM 2 
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^±{r±,nt ot ) . (A.2) 



where f2 r f = [i B B^j2^2 is the coupling strength (taking the Lande factor g F = —1/2). 
Vb = V g + V»b is the sum of the transverse potential V g = \muJ\r\ and the longitudinal 
confinement V\\ t B = \muj 2 z 2 . As discussed in section l2~2l we write the total wavefunction 
in the separate form 

<f>h(z,t) 

with the normalization condition J \ip±\ 2 dr± = 1. The global normalization is here 
J \ij)\ 2 dr±dz = N such that |0b| 2 = ^b and |0l| 2 = n L correspond respectively to the 
BEC and atom laser linear density. n tot = n^ + n^ is the total linear density. Following 
the procedure of reference [24] (i.e. minimizing the action associated to the functional 
energy (ij:\'H\^)), one obtains first the transverse stationary GP equation 
/ h 2 \ 

[ - 7^ A ± + V S + gntat\i>±\ 2 J 4>± = /ieffKot) i/>± , (A.3) 

where /x e ff is the effective local chemical potential. From equation ( 1A.3[) . ip± depends 
implicitly on z via the total linear density, as written in (1A.2j) . This property results 
from the hypothesis on g being independent of the magnetic substates. 

In the so-called "ID mean-field" regime (an tot < 1), the weak interactions do not 
perturb ip±, which sticks to the fundamental gaussian mode and yU e g ~ hui±(l + 2an tot ). 
In the opposite strong density regime (3D TF regime defined by an tot 3> 1), ip± is 
well approximated by an inverted parabola and /i e g ~ 2huj_y/an tot . Altogether the 
phenomenological expression 

/i e ff = fkJ±y/l + Aan tot (A.4) 

is very accurate for any interaction strength [38]. Using equation f ]A.4j) . one finally 
derives the quasi-lD coupled GP equations for the longitudinal dynamic: 

^"HT = f"7T~~ ^-^rf+V|| ) B+/i e ff)</'B--?^0L (A.5) 
ot \2m J 2 

The numerical simulations presented in figure [7] are obtained from this set of equations. 
Appendix A.2. Approximated form 

The above set of equations can be simplified assuming i) the adiabatic following of the 
BEC wave function and ii) no intra-laser interactions (an^ = 0). In the TF regime, 
equation ( 1A.5I) reduces indeed to 

^bec = ^mu 2 z z 2 + hwj_^/l + 4an B (z) , (A. 7) 

and the BEC linear density writes 

I r/ /iBEc \ 2 /, _ z^ 2 
XhujJ \ RV 



n B[ z) = ^ 



- 1 



e(i-iS-). (A.s) 
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The BEC edge is located at R z = (2/2 B Ec/ mw D 1/ ' 2 ~ Rz ~ ^_i_/2/z B ec> with /ibec — 
yU-BEC — hu)±- Except around this position (where the density is extremely weak), it is 
very well approximated by Ub{z) ~ 15iVo/16i? 2 (l — z 2 / R 2 ,) 2 (z <C R z ), leading to the 
usual expression /xbec — {15aNQh 2 ujj_uj z ) 2 ' 5 'm 1 ' 5 /2. Here the BEC depletion has been 
neglected but the time dependence can be included using Ns(t) — N — 2Al(£) (where 
TVl is the atom number for one side) instead of N in the preceding expressions. 
From equation flA.41) . the effective chemical potential reads finally 



z 2 \ „ / : 



fi eS (z) =Hu ± + j2 BEC (1 - O (1 - 4^ J . (A.9) 

It acts as an effective ID potential for the atom laser (V\\ = fi e s), whose dynamics is 
described by the Schrodinger equation 

For dilute atom laser beams, the intra-laser interactions can be added in a 
perturbative way in the expression ( 1A.4I) . In the guide, it corresponds to the ID mean- 



field regime (an L <C 1). There fi e s ~ hu±(l + 2anj J ) and the propagation of 0l along 
the guide (with the eventual presence of an obstacle Kbst) is governed by the nonlinear 
equation 

ih-^ = ^ - — — + V ohst + 2frw±an L J L . (A.ll) 

Appendix B. Stationary longitudinal atom laser wave function </>l,e 

We calculate here the atom laser wave function l ,e at the energy E = /Ubec — Mui^. 
Those are the stationary solutions of the ID Schrodinger equation in the effective 
potential V\\(z) [see equation flA.lOj) ]. In dimensionless unity (length and energy units 



being respectively a z and hu z /2), it writes 

1 [z) = 0, (B.l) 



~ - z 2 e(l-\z\/R z ) + 5 



where 5 = 25u r f/u z is the frequency detuning and R z = R z /a z the corrected TF radius. 
In the BEC zone (z < R z ), the potential is a inverted harmonic oscillator and the odd 
and even stationary solutions are [391 HO] 

Mz) = ze^ 2/2 F( 3 -^\ 3 -\iz 2 ), (B.2) 

c (,-) = e-- 2 / 2 F(i^|i|iz 2 ), (B.3) 

where F(a|b|x) are the confluent hypergeometric functions |41j . The BEC wave function 
being symmetric, the odd contribution vanishes in the overlap integral of the Fermi 
golden rule [equation (j3J)]. We thus only consider the even function (0 [z < R z ) = 
Acj) e (z), with A a prefactor to be determined below) and restrict the analysis to z > 0. 
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Figure Bl. Longitudinal wave function 4>(z) for different normalized rf detunings 
(S = —5, 0, 20 and 50). The numerical parameters (/2bec = ^P-bec/^z = 185, 
corresponding to R z = 13.6) are those of section [3~4l zb = S 1 ^ 2 is the classical turning 
point position in normalized unit. 



Outside the BEC, the wave function <p is connected to a stationary plane wave. 
Introducing a quantization box of size L, it takes the form 

1 



z >R Z 



sin(fcz + if) 



(B-4) 



with the dimensionless momentum k = E 1 / 2 associated to the atom laser kinetic energy 
in the guide. From the continuity of <fi and its derivative at the BEC edge (z = R z ), the 
two free parameters A and (p write 



A 



VI F(a|l/2|iR z 2 )' 



— kR z + arctan 



ik/R z 



_ a F(a+l|3/2|iB|) 
1 ^ U F(a|l/2|iR2) 



(B.5) 



(B.6) 



a/6)F(a+l|b + l|x) [H]. Note 
, apart around the classical 



Here a = (1 — i5)/A and we use the relation ^(a\b\x) = 
that a WKB calculation will lead to a similar result 
turning point z E . 

From the knowledge of 0l,e (see figure |BT|) and 0b [equation flA.8[) ] . we calculate the 
overlap integral in the outcoupling rate expression ([3]). Here the density of states for a 
one-dimensional box of size L is rj(E) = L(2m/E) 1 / 2 /2%h (the presence of the mean-field 
potential in the overlap region with the BEC has a negligible effect for L — > oo). 
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